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Abstract — We  investigate  the  total  resistive  losses  incurred 
in  returning  a  power  network  of  identical  generators  to  a 
synchronous  state  following  a  transient  stability  event  or  in 
maintaining  this  state  in  the  presence  of  persistent  stochastic 
disturbances.  We  formulate  this  cost  as  the  input-output  H2 
norm  of  a  linear  dynamical  system  with  distributed  distur¬ 
bances.  We  derive  an  expression  for  the  total  resistive  losses 
that  scales  with  the  size  of  the  network  as  well  as  properties 
of  the  generators  and  power  lines,  but  is  independent  of  the 
network  topology.  This  topologically  invariant  scaling  of  what 
we  term  the  price  of  synchrony  is  in  contrast  to  typical  power 
system  stability  notions  like  rate  of  convergence  or  the  region 
of  attraction  for  rotor-angle  stability.  Our  result  indicates  that 
highly  connected  power  networks,  whilst  desirable  for  higher 
phase  synchrony,  do  not  offer  an  advantage  in  terms  of  the 
total  resistive  power  losses  needed  to  achieve  this  synchrony. 
Furthermore,  if  power  flow  is  the  mechanism  used  to  achieve 
synchrony  in  highly-distributed-generation  networks,  the  cost 
increases  unboundedly  with  the  number  of  generators. 

I.  Introduction 

Many  factors  such  as  increased  demand,  renewable  energy 
mandates  [1]  and  further  deregulation  of  the  electric  power 
industry  [2],  [3],  [4]  are  driving  changes  to  the  electric  power 
grid.  The  new  grid  will  have  to  deal  with  higher  levels  of 
uncertainty  from  renewable  energy  sources  and  changing 
load  patterns  as  well  as  increasingly  distributed  energy 
generation.  These  changes  are  likely  to  make  it  more  prone 
to  stability  issues.  In  particular,  they  have  the  potential  to 
create  problems  associated  with  rotor-angle  stability,  which 
is  the  ability  of  the  power  grid  to  recover  synchrony  after  a 
disturbance  [5],  Synchrony  refers  to  the  condition  when  both 
the  frequency  and  phase  of  all  generators  within  a  particular 
power  network  are  aligned.  Loss  of  synchrony  can  lead  to 
load  shedding  or  even  black-outs. 

A  special  case  of  rotor-angle  stability  is  the  so-called 
transient  stability  problem,  which  is  associated  with  large 
angle  disturbances  due  to  events  such  as  generator,  power 
line  or  other  component  failures  or  other  abrupt  changes  that 
can  be  caused,  for  example,  by  intermittent  renewable  energy 
sources.  There  is  a  large  body  of  power  system  transient 
stability  literature,  see  e.g.  [6],  [7],  Most  of  this  literature 
focuses  on  the  existence  of  Lyapunov  like  energy  functions 
[8],  [9]  and  their  use  in  determining  a  region  of  attraction 
type  criteria  for  a  particular  synchronous  state  or  set  of  states, 
see  e.g.  [10],  [11], 
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A  recent  research  trend  has  been  to  draw  connections 
between  problems  in  distributed  system  control  and  power 
network  stability.  This  literature  is  vast,  but  we  note  in 
particular  a  series  of  works  that  use  a  set  of  coupled  Ku- 
ramoto  oscillators  to  represent  the  power  network  [12],  [13]. 
This  non-uniform  Kuramoto  oscillator  modeling  framework 
uses  a  first  order  approximation  of  the  network  reduced 
classical  power  system  model  to  provide  network  parameter 
dependent  analytical  conditions  for  frequency  and  phase 
synchronization  [13].  Similar  first  order  models  have  been 
employed  to  investigate  how  power  flow  scheduling  and 
adding  more  power  lines  (i.e.,  increasing  graph  connectivity) 
affects  the  rate  of  convergence  in  a  power  network  [14]. 

In  the  present  paper,  we  examine  the  connections  between 
distributed  control  systems  and  power  network  stability  in  a 
different  context.  We  do  not  study  stability,  but  rather  assume 
that  the  network  will  return  to  a  synchronized  state  after 
disturbances.  Instead  we  focus  on  the  cost  of  keeping  the 
network  in  synchrony,  i.e.  how  much  real  power  is  required  to 
drive  the  system  to  a  stable,  synchronous  operating  condition. 
Lack  of  synchrony  leads  to  circulating  currents  [15]  passing 
between  generators  whose  angles  are  out  of  phase.  This  flow 
of  current  leads  to  resistive  power  losses  over  the  power 
lines  due  to  their  non-zero  line  resistances.  This  loss  is 
generally  considered  relatively  small  compared  to  the  total 
real  power  flow  in  a  typical  power  network.  It  is  however 
unclear  whether  these  losses  will  be  small  in  power  grids 
of  the  future,  which  are  expected  to  have  highly  distributed 
generation,  and  consequently  many  more  (though  typically 
smaller  power  capacity)  generators  than  today’s  grid. 

The  problem  we  analyze  is  that  of  a  large  network  of  many 
identical  generators.  We  consider  several  scenarios  such 
as  the  power  network  encountering  disturbance  (transient 
stability)  events,  or  being  subjected  to  persistent  stochastic 
disturbances  where  the  system  is  continuously  correcting  for 
these  disturbances.  In  both  of  these  scenarios,  we  quantify  the 
total  power  lost  due  to  non-zero  line  resistances.  Our  main 
result,  in  equation  (13),  shows  that  this  power  loss  scales 
with  the  product  of  the  network  size  and  the  ratio  of  line 
resistances  to  their  reactances  (which  is  assumed  to  be  the 
same  for  all  links).  The  latter  quantity  is  normally  assumed 
to  be  rather  small.  However,  our  result  shows  that  even 
though  that  ratio  is  small,  when  the  number  of  generators 
becomes  very  large,  then  so  will  the  total  resistive  losses. 
Furthermore,  and  perhaps  more  surprisingly,  these  losses  are 
independent  of  the  network  topology,  i.e.  highly  connected 
networks  (which  would  have  highly  coherent,  or  close  to 
identical,  phases)  and  loosely  connected  networks  (which 
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would  have  a  higher  degree  of  phase  incoherence)  incur  the 
same  resistive  power  losses  in  recovering  synchrony. 

Our  results  indicate  that  the  cost  of  maintaining  synchrony 
using  power  flows  is  essentially  a  function  of  the  number  of 
generators  in  the  network  and  not  of  its  topology.  We  should 
point  out  that  this  conclusion  is  not  inconsistent  with  other 
results  on  power  system  stability  and  performance  measures. 
For  example,  interaction  strength  and  network  topology  play 
important  roles  in  determining  whether  a  system  can  syn¬ 
chronize  [8],  [12],  [13],  [16]  and  the  rate  of  convergence  or 
damping  of  a  power  system  is  directly  related  to  the  network 
connectivity  [14],  The  numerical  examples  we  present  in 
this  paper  confirm  that  while  losses  are  independent  of 
the  network  connectivity,  a  highly  connected  network  will 
stabilize  more  quickly  and  with  less  oscillatory  behavior. 

One  measure  that  quantifies  the  lack  of  synchrony  between 
network  elements  is  that  of  “coherence”  [16],  As  expected, 
highly  connected  networks  are  more  coherent  than  loosely 
connected  ones.  We  refer  the  reader  to  [16]  for  asymptotic 
scalings  of  lattice-type  networks  and  to  [17]  for  fractal 
networks,  and  note  that  these  scalings  are  somewhat  different 
than  those  that  quantify  network  damping.  One  intuitive 
explanation  for  the  cost  of  synchrony  being  independent 
of  network  topology  is  as  follows.  Comparing  a  highly 
connected  network  with  a  loosely  connected  one,  we  expect 
the  former  to  have  much  more  phase  coherence  than  the 
latter.  Consequently  the  power  flows  per  link  in  a  highly 
connected  network  are  relatively  small,  but  there  are  many 
more  links  than  in  the  loosely  connected  network.  Thus  in  the 
aggregate,  the  total  power  losses  of  the  two  networks  are  the 
same.  One  should  keep  in  mind  however,  that  a  less  coherent 
network  is  more  likely  to  have  transient  stability  problems, 
exit  the  region  of  attraction,  etc.  The  issues  of  stability  and 
the  cost  of  synchrony  are  two  different  issues. 

The  remainder  of  this  paper  is  organized  as  follows. 
Section  II  introduces  the  classical  power  system  model  and 
defines  the  notation.  In  Section  III  we  quantify  total  resistive 
line  losses  as  an  input-output  H 2  norm  of  a  linear  dynamical 
system.  This  framework  is  then  used  to  derive  an  algebraic 
expression  for  the  resistive  line  losses  in  terms  of  the  param¬ 
eters  of  the  admittance  matrix,  generator  damping  and  the 
size  of  the  network.  Various  interpretations  of  the  II2  norm 
are  then  presented  together  with  the  corresponding  operating 
scenarios  for  a  power  network  .  The  results  of  the  analysis  are 
demonstrated  for  two  different  network  topologies  in  Section 
IV.  The  paper  concludes  and  suggests  directions  for  future 
work  in  Section  V. 

II.  Problem  Formulation 

Consider  a  network  of  N  buses  (nodes)  and  £  edges.  At 
each  node  i  =  1. ....  A  there  is  a  generator  6’,,  with  inertia 
constant  Mj,  damping  6t,  voltage  magnitude  V-L  and  angle 
Using  the  classic  machine  model,  which  assumes  a  fixed 
voltage  behind  a  reactance,  the  dynamics  of  the  ith  generator 
is  given  by  [18] 

Midi  +  pdi  =  Pm,i  -  Pe,i  v  i  =  1,  2 ...n.  (1) 


Gi  G,  G3  Gn.,  G„ 

ji;  A]  [KA]  {r„s,}  feiilfer  {KA} 

u 

£  Pi  p,  p.-,  p. 

Fig.  1:  A  depiction  of  a  linear  network  of  n  generator  nodes 
(buses)  Gi  for  i  =  1 . . . ,  n. 

Here,  Pe  i  is  the  electrical  power  flowing  out  of  the  ith 
generator,  and  PTOi;  is  the  mechanical  power  input  from  the 
corresponding  turbine. 

Define  the  admittance  over  edge  £tj  connecting  nodes  i 
and  j  as  ytJ  =  gij  —  \bl:j,  where  gtj  and  bi:j  are  respectively 
the  conductance  and  susceptance  of  the  line  defined  by  the 
edge  £ij.  Then,  the  electrical  power  injection  at  node  i  is 

Pe,i  =  gu\Vj\2  +  y^,gik\Vj\  \Vk\  cos (Oj  -  Ok) 

kr^i 

+  J2bik\yt\\Vk\sm(0i^0k),  (2) 

kr^i 

where  k  ~  i  denotes  an  edge  £lk  and  gu  is  the  self 
admittance  of  the  generator.  In  what  follows,  we  use  a 
simplified  model  in  which  all  generators  are  assumed  to  have 
the  same  self  admittance,  which  we  denote  g  =  glt  for  all 
i  =  1, . . . ,  n.  Finally,  applying  the  standard  linear  power  flow 
assumptions  used  in  transient  stability  analysis  (conductances 
are  negligible,  angle  differences  are  small  and  voltages  are 
constant  with  unit  magnitude)  to  (2)  yields 

Pe,i  »  53  bik  [Oi  -  0k\ .  (3) 

fceA 

Substituting  this  into  (1)  leads  to 

MA  +  pOi  »  -  53  btk  [Oi  -  0k]  +  Pm,j.  (4) 
fee  A 

In  order  to  simplify  the  notation  we  define  the  entries  of 
the  admittance  matrix  Y  G  Cnx™  as 

{5 +  53  (/iP  if*  =  J. 

- Vij ,  if  i  ±  j  and  j  ~  z,  (i.e.,  %  G  £), 

0  otherwise. 

Then,  we  partition  this  admittance  matrix  into  the  real 
(resistive)  and  imaginary  (reactive)  parts  and  define 

y  =  Rem  +  jim m  -  m+^)  +  j  lb. 

By  this  construction,  Lq  and  L  n  retain  the  symmetry  of  Y, 
and  they  have  as  a  common  eigenvector  the  vector  1  with 
components  all  equal  to  1,  i.e. 

Lb  1  =  Lq  1  =  0. 

It  is  a  well-known  result  that  if  the  graphs  underlying  the 
system  represented  by  and  Lq,  which  are  analogous 


to  weighted  graph  Laplacians,  are  connected  (i.e.  any  two 
nodes  are  connected  by  a  path  of  edges),  then  the  remaining 
eigenvalues  of  Lb  and  LG  are  all  positive. 

An  assumption  we  will  invoke  later  is  that  all  the  eigen¬ 
vectors  of  Lb  and  Lq  are  shared  (u  is  an  eigenvector  of  L /> 
iff  it  is  an  eigenvector  of  Lb)-  An  important  consequence  of 
this  assumption  is  that  the  two  Laplacians  are  simultaneously 
diagonalizable  by  the  same  orthogonal  transformation  U,  i.e. 

U*LbU  =  A  s,  U*LgU  =  A  G,  (5) 

where  A b  and  Aq  are  diagonal  matrices  with  the  respective 
eigenvalues  of  Lb  and  Lq  as  diagonal  entries.  One  setting 
in  which  this  assumption  holds  is  when  the  ratios  of  each 
connection’s  conductance  to  susceptance  (equivalently  the 
ratios  of  resistances  to  reactances)  are  all  equal,  i.e. 

9ik  '  ik 

j—  =  —  =  a, 

®ik  -kik 

which  is  independent  of  the  link  indices  ( i ,  k).  It  then  follows 
that 

Lq  =  aLBl  (6) 


which  implies  that  LG  and  Lb  have  the  same  eigenvectors. 
Finally,  we  rewrite  equation  (4)  in  state  space  form 


d 

' m ' 

0  I 

9(t) 

dt 

w(i)_ 

-Lb  -£I_ 

w(f)_ 

=:  Aij)  +  B  w, 


where  we  have  assumed  that  the  constant  M*  =  1  for  all  i 
and  that  mechanical  power  input  Pm y  is  a  constant  that  can 
be  captured  as  a  part  of  the  input  w. 

Remark  1:  Based  on  the  classical  machine  model  F>rn  t  is 
essentially  the  real  power  input  at  the  bus.  Thus  one  can 
consider  Pmy  as  a  system  input.  Alternatively,  the  absence 
of  Pm,i  can  be  justified  as  follows.  For  a  system  with  no 
loads  the  equilibrium  point  of  the  system  is  at  Pey  =  Pmy. 
If  (1)  defines  a  linear  system  (as  assumed  here)  a  coordinate 
transform  can  be  applied  to  obtain  a  new  system  with  shifted 
angles  9s  corresponding  to  the  transformed  system  with  the 
equilibrium  point  Pm,i  =  0.  By  abuse  of  notation  we  denote 
the  shifted  coordinates  as  9  in  the  state  space  representation 

(7). 

Remark  2:  The  network  model  (7)  assumes  a  system  with 
no  loads.  Such  a  model  may  arise  from  network  reduction 
leading  to  a  system  of  internal  generator  buses  with  loads 
approximated  as  impedances  that  have  been  absorbed  into  the 
network.  The  model  could  easily  be  extended  to  include  static 
loads  through  a  system  of  differential  algebraic  equations 
similar  to  those  derived  in  [14],  [19]. 


III.  System  Performance 

Several  performance  metrics  can  be  used  to  quantify  the 
relative  stability  for  the  system  (7).  A  few  common  metrics 
are  the  ability  to  synchronize,  the  degree  of  synchronization 
that  is  achievable,  the  time  to  synchronize,  or  the  control 
effort  required  to  obtain  the  desired  system  state.  In  the 
distributed  systems  setting  it  is  also  common  to  evaluate 
these  metrics  with  respect  to  various  control  strategies.  For 


example,  in  evaluating  the  performance  of  a  distributed 
system  with  local  versus  global  control  strategies.  In  this 
section  we  focus  on  the  control  actuation  required  to  drive 
a  system  to  the  synchronous  state.  This  control  effort  is 
measured  through  the  real  (resistive)  power  loss  over  each 
line.  These  losses  are  associated  with  circulating  currents 
that  arise  from  the  angle  differences  between  generators,  i.e. 
disturbances  [15], 

The  power  flow  over  an  edge  £,3  is 

Pij  +  Pji  =  Vi  (Vi  -  Vj)*  Vij  +  Vj  (Vj  -  Vi)*  y„. 

where  *  denotes  the  complex  conjugate.  The  resistive  power 
loss  over  £tJ  can  therefore  be  defined  as 

P%ss-.=  gik\Vi-Vk\2.  (8) 

Using  a  small  angle  approximation  and  standard  trigonomet¬ 
ric  identities  this  can  be  approximated  as 

P!r  =9ik\0i-9k\2.  (9) 

We  are  interested  in  the  sum  total  of  the  resistive  losses  over 
all  links  in  the  network,  which  is  given  by 

Pioss  =  Y,9ik\0i-ek\2  -  (io) 

This  last  quantity  can  be  expressed  in  vector  form  as 

P;oss  =  y*y 


where  the  vector  signal  y  is  an  output  of  the  linearized  swing 
equations  (7) 


y  =  Cty 


with  Ci Ci  :=  LG-  (11) 


A  simple  choice  for  C\  would  be  to  take  C \  =  Lq  (possible 
since  LG  is  positive  semi-definite),  which  is  what  we  assume 
from  now  on. 

We  now  calculate  the  H2  norm  from  disturbance  w  to 
the  performance  output  y  of  the  system  (7)  and  (11).  The 
square  of  the  H2  norm  has  several  standard  interpretations 
including  (a)  The  variance  of  the  output  y  when  the  input 
w  is  a  unit  variance  white  stochastic  process,  (b)  The  total 
time  integral  of  the  variance  of  y  when  the  initial  condition 
is  a  random  variable  with  correlation  matrix  BB*,  and  (c) 
The  total  sum  of  time  integrals  of  output  response  powers 
given  an  impulse  as  a  disturbance  input  at  each  generator. 

We  first  perform  the  H2  norm  calculation  and  derive  a 
formula  in  terms  of  the  system’s  parameters,  and  then  inves¬ 
tigate  the  implications  of  the  three  different  interpretations 
for  this  particular  system  of  swing  equations. 


H2  Norm  Calculation 


For  ease  of  reference  we  rewrite  the  system  equations  (7) 
and  (11)  here 
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We  will  denote  the  input-output  mapping  of  this  system 
by  H.  This  system  has  all  eigenvalues  strictly  in  the  left 
half  of  the  complex  plane  with  the  exception  of  a  single 
zero  eigenvalue  of  Lb-  It  is,  however  easy  to  show  that 
this  unstable  mode  is  unobservable  from  the  output  y  due 
to  the  fact  that  Lq  =  Cf  Cj  shares  this  eigenvalue  and 
its  corresponding  eigenvector  (see  Appendix  for  the  full 
argument).  It  then  follows  that  the  input-output  transfer 
function  from  w  to  y  is  indeed  stable  and  has  finite  H 2 
norm. 

The  easiest  method  to  calculate  the  H2  norm  of  the  (12) 
is  using  a  spectral  decomposition  of  Lb-  Consider  the  state 
transformation 
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where  U  is  the  matrix  in  (5)  diagonalizing  Lb  and  Lq-  The 
system  in  the  new  state  variables  becomes 


I 

~P I 


Since  multiplying  by  orthogonal  matrices  does  not  change 
the  H 2  norm,  we  can  multiply  the  input  by  U  and  and  the 
output  by  U*  (i.e.  define  w  =  17*  w  and  y  =  U*y)  to  obtain 
an  equivalent  system  (that  has  the  same  H 2  norm) 
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We  will  denote  the  input-output  mapping  of  this  system 
by  H.  Since  A b  and  Aq  are  diagonal,  this  represents  N 
decoupled  systems 
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where  n  =  1 , ,N  are  the  indices  of  eigenvalues  Af  and 
Af ,  which  correspond  to  Lb  and  Lq  respectively.  Denote 
the  input-output  mapping  of  each  decoupled  subsystem  by 
.  We  can  then  write 


H  =  diag  TTjv)  . 


The  square  of  the  H2  norm  of  the  system  (12)  is  thus  the 
sum  of  the  squares  of  the  H2  norms  of  all  the  decoupled 
subsystems 


N 

ml  =  wmi  = 

n= 1 

The  H 2  norm  of  each  of  the  individual  subsystems  can  now 
be  easily  calculated  as  follows.  For  n  =  1,  the  corresponding 
eigenvalues  are  Af  =  Af  =  0,  and  we  have  a  completely 


unobservable  system,  thus  ||ffi||2  =  0.  For  rt  ^  1,  we  solve 
the  Lyapunov  equation  for  the  observability  Grammians  Xn 
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where  for  simplicity  of  notation  we  have  dropped  the  sub¬ 
script  n  from  the  components  of  of  the  Grammian  Xn.  This 
matrix  equation  corresponds  three  equations,  of  which  only 
the  following  two  are  relevant 

-Af  X*  -  X0Af  =  -  Af  ^Re(X0)  =  ^ 

Xq  +  X*  -  2(3X22  =0  ^X22  =  ^Re(X0). 


Finally,  since  the  B  matrix  of  each  subsystem  is  [0  1]T,  the 
H2  norm  (squared)  of  each  subsystem  is  just  X22,  and  we 
therefore  conclude  that 


II 4.  II 
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The  total  H 2  norm  (squared)  of  the  overall  system  (12)  is 
thus 
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In  summary,  we  conclude  that  the  total  resistive  losses  in 
this  network  are 
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(13) 


where  0  is  a  generator’s  self  damping,  -  is  the  ratio  of  a 
line’s  resistance  to  its  reactance  (assumed  equal  for  all  lines), 
and  N  is  the  number  of  generators  in  the  network.  Note  that 
this  expression  is  independent  of  the  network  topology  but 
depends  on  the  number  of  generators  in  the  network. 


H 2  Norm  Interpretations  for  Swing  Dynamics 

In  the  previous  subsection  we  calculated  the  H2  norm  of 
the  linearized  swing  dynamics  (7),  together  with  the  output 
equation  (11)  based  on  a  disturbance  input  (forcing)  w.  In 
our  formulation  the  square  of  the  Euclidean  norm  y*y  of  the 
output  vector  was  defined  to  be  equal  to  the  total  real  power 
dissipated  due  to  resistive  losses  in  the  network  connections 
as  the  system  synchronized.  This  synchronized  state  (the 
equilibrium)  is  maintained  by  circulating  currents  leading  to 
power  flowing  back  and  forth  between  generator  nodes. 

It  is  well  known  that  the  H 2  norm  has  at  least  three 
different  interpretations.  We  recall  these  interpretations  here 
in  order  to  give  three  different  physical  scenarios  in  which 


(a)  5  Node  Line  Graph  (b)  Fully  Connected  5  Node  Graph 


(c)  20  Node  Line  Graph  (d)  Fully  Connected  20  Node  Graph 


Fig.  2:  (a),(b)  Simulation  of  a  5  bus  system  of  identical  oscillators  with  a  random  initial  velocity  disturbance,  and  zero 
initial  phase  disturbance.  (c),(d)  Simulation  of  a  system  with  20  identical  oscillators  under  the  same  conditions  as  (a),(b). 
Figures  (a),(c)  are  for  linearly  connected  networks  (with  line  graphs  as  in  Figure  1),  while  Figures  (b),(d)  are  for  fully 
connected  networks  (complete  graphs).  Note  how  linear  networks  are  much  less  coherent  than  complete  networks  (Figures 
(a),(c)  versus  (b),(d)  respectively).  Total  resistive  losses  are  however  the  same  for  the  5  node  networks  (a)  and  (b),  and  for 
the  20  node  networks, (c)  and  (d),  respectively,  reflecting  that  resistive  line  losses  are  a  function  of  the  network’s  size  rather 
than  its  topology. 


equation  (13)  quantifies  the  resistive  losses  due  to  the  syn¬ 
chrony  requirement.  Denote  by  H  the  following  Linear  Time 
Invariant  (LTI)  system 

ip(t)  =  Aip(t)  +  Bw(t) 

V(t)  =  Cip(t), 

The  H2  norm  ||fT||2  of  the  system  has  the  following  three 
different  interpretations. 

(a)  Response  to  a  white  stochastic  input.  When  the  input 
w  is  a  white  second  order  process  with  unit  covariance 
(i.e.  E  {w(t)w*(t)}  =  I),  the  H 2  norm  (squared)  of  the 
system  is  the  steady-state  total  variance  of  all  the  output 
components,  i.e. 

\\H\\l  =  lim  E  {y* (t)y(t)}  . 

£—>•00 

For  the  swing  dynamics  (12)  the  disturbance  vector 
can  be  thought  of  as  persistent  stochastic  disturbances 
(forcing)  at  each  generator.  These  disturbances,  which 


are  uncorrelated  across  generators,  can  be  due  to  un¬ 
certainties  in  local  generator  conditions,  sudden  changes 
in  load,  or  fault  events.  The  variance  of  the  output  is 
exactly  the  expectation  of  the  total  power  loss  due  to 
line  resistances. 

(b)  Response  to  a  random  initial  condition.  When  the  input 
is  turned  off,  but  the  initial  condition  is  a  random  variable 
x0  with  correlation  E  {x0x*0}  =  BB *,  then  the  H 2  norm 
is  the  time  integral 

pOO 

ml  =  /  E{y*(t)y(t)}  dt 

J  o 

of  the  resulting  response  y. 

The  interpretation  for  (12)  is  as  follows.  Since  BB*  = 

Q  ^  ,  the  corresponding  initial  condition  corresponds 

to  each  generator  having  a  random  initial  velocity  per¬ 
turbation  that  is  uncorrelated  across  generators  and  zero 
initial  phase  perturbation.  Then  \\H\\l  quantifies  the  total 


(over  all  time)  expected  resistive  power  losses  due  to  the 
system  returning  to  a  synchronized  state. 

(c)  Sum  of  responses  to  impulses  at  all  inputs.  Let  ej  refer  to 
the  vector  with  all  components  zero  expect  for  1  in  the 
ith  component.  Consider  N  experiments  where  in  each, 
the  system  is  fed  an  impulse  at  the  ith  input  channel, 
i.e.  w;(t)  =  ej(5(t).  Denote  the  corresponding  output 
by  yi.  The  If2  norm  (squared)  is  then  the  sum  total  of 
the  L 2  norms  of  these  outputs,  i.e. 

N  pOO 

\\h\\2  =  &*(%»(*)  dt. 

i= i Jo 

A  stochastic  version  of  this  scenario  corresponds  to 
a  system  where  the  inputs  w;  can  occur  with  equal 
probability.  Under  this  assumption  the  \\H\\\  becomes 
the  expected  total  power  loss  given  these  inputs. 

The  interpretation  for  (12)  is  when  each  of  the  generators 
are  subject  to  impulse  force  disturbances  (since  enters 
into  the  momentum  equation  of  each  generator),  and 
||TZj| 2  is  thus  the  total  power  loss  over  all  time  under 
such  a  scenario. 

IV.  Numerical  Examples 

Consider  two  networks  of  identical  generators,  one  whose 
underlying  graph  is  a  line,  as  in  Figure  1,  and  one  with 
a  fully  connected  graph.  In  this  section,  we  compare  the 
behavior  of  systems  with  these  two  network  topologies  as 
the  network  size  varies.  All  simulations  use  the  following 
parameter  values  [20]:  M  =  /?  =  with  a  frequency 

f  =  60  Hz.  The  admittance  between  connected  generators 
is  set  to  i/ij  =  0.2  +  jl.5  and  we  assume  gti  =:  g  =  0  for 
all  i  =  1, . . . ,  n. 

Figures  2a  and  2b  respectively  show  the  state  trajecto¬ 
ries  of  a  5  node  system  with  a  line  graph  and  one  with 
a  fully  connected  graph  given  identical  initial  conditions. 
The  initial  conditions  for  each  Wj(0)  were  drawn  from  a 
uniform  distribution  [—1,1],  and  0(0)  was  set  to  zero.  This 
corresponds  to  the  H2  norm  interpretation  (b)  described 
in  the  previous  section.  For  both  the  line  graph  and  the 
fully  connected  network  the  total  losses  are  the  same,  with 
p/oss  _  i  2885,  as  predicted  by  the  main  result  in  (13). 
However,  the  transient  behavior  shows  that  the  system  with 
the  fully  connected  graph  has  more  “coherent”  phases  (i.e. 
they  stay  closer  together).  This  is  most  clearly  visible  in 
the  top  panels  of  the  plots  that  depict  the  state  0.  Here,  the 
small  oscillations  for  the  line  graph  system  continue  through 
the  10  second  interval  shown  in  Figure  2a  but  have  almost 
completely  died  out  for  the  fully  connected  of  Figure  2b. 

In  order  to  evaluate  the  effect  of  increasing  the  network 
size  we  ran  a  similar  test  with  two  20  bus  systems  with  the 
same  topological  structure  and  initial  conditions  (uii(0)  £ 
[—1,1]  drawn  from  a  uniform  distribution  for  each  i  and 
0(0)  =  0).  In  this  case,  the  losses  increased  to  P,oss  = 
5.9256  but  remained  equal  in  both  networks.  The  faster 
convergence  to  the  synchronized  state  in  fully  connected 
system  (compared  to  the  line  graph)  is  much  more  evident 
in  the  larger  network,  as  shown  in  figures  2c  and  2d. 


The  simulation  results  indicate  that  a  fully  connected 
graph  is  much  more  coherent  than  a  line  graph.  The  resem¬ 
blance  of  this  result  to  those  about  coherence  in  vehicular 
formations  (platoons)  and  similar  consensus-like  network 
algorithms  [16]  is  striking.  With  coherence  as  a  performance 
measure,  the  line  graph  is  the  worst  such  topology  (amongst 
connected  graphs),  while  the  fully  connected  topology  is  the 
best.  However  this  additional  coherence,  while  desirable  for 
other  reasons,  will  not  result  in  lower  resistive  losses. 


V.  Conclusions 

We  have  considered  a  power  network  model  with  dis¬ 
tributed  disturbances,  and  quantified  the  total  resistive  power 
losses  incurred  due  to  the  current  flows  needed  to  maintain 
phase  synchrony.  We  have  shown  that  these  losses  are  in¬ 
dependent  of  the  network’s  topology,  but  scale  unboundedly 
with  the  number  of  generators.  There  are  interesting  impli¬ 
cations  for  the  design  of  future  highly-distributed-generation 
networks,  which  have  potentially  orders  of  magnitude  more 
generator  nodes  that  today’s  networks.  While  a  highly  phase 
coherent  (and  therefore  highly  connected)  network  is  desir¬ 
able  for  many  reasons,  the  cost  of  maintaining  this  coherence 
depends  only  on  the  number  of  generators  and  not  the 
network’s  connectivity.  Since  this  cost  grows  unboundedly 
with  the  number  of  generator  nodes,  then  the  current  scheme 
of  using  power  flows  as  the  synchronization  mechanism  may 
not  be  scalable  to  future  networks.  This  is  perhaps  a  further 
argument  for  the  use  of  other  control  mechanisms,  such  as 
communication  links,  for  phase  synchronization. 


Appendix 


Proof  that  (12)  is  input-output  stable 


It  is  well  known  that  for  any  pair  (C,  A),  the  observability 
of  (C.  A)  and  (C*C.  A )  are  equivalent.  The  only  unstable 
mode  of  the  A-matrix  in  (12)  is  at  zero  with  corresponding 


eigenvector 
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.  This  mode  is  unobservable  from  the  output 


y  since  by  the  PBH  test 
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Note  that  this  is  a  consequence  of  Lq  and  Lb  having  the 
common  eigenvector  1  with  eigenvalue  0.  The  unobservabil¬ 
ity  of  the  only  unstable  mode  implies  then  that  the  system 
(12)  is  input-output  stable. 
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